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DEPENDENCE OF THE VISCOSITY OF BASALT MELTS ON THE 
CHEMICAL COMPOSITION OF THE INITIAL MINERAL MATERIAL 
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An approximate relation for calculating the viscosity of fused basaltic glass according to the chemical compo¬ 
sition of the initial mineral raw material is proposed. This relation was obtained as a result of a multifactor re¬ 
gression analysis of a set including 27 rocks. The experimental values of the viscosity of melts and x-ray fluo¬ 
rescence data on the content of oxides in the raw material were used in the modeling. The relation makes it 
possible to predict to a high degree of accuracy the viscosity of melts at a prescribed temperature. 

Key words: rocks, mineral raw material, chemical composition, acidity modulus, melt, viscosity, regression 
analysis. 


One of the most important physical-chemical properties 
of aluminosilicate glass that determine the suitability of the 
glass for manufacturing different types of fibers is the visco¬ 
sity r|, which is usually compared with the values of the 
acid-base indicators, first and foremost, the acidity modulus 
M a calculated as the ratio of the mass content (%) of four 
melt-forming oxides: 

w Si0 2 +Al 2 0 3 

M n = ---. 

a CaO+MgO 

The larger the value of M a of the initial raw material, the 
higher the viscosity of its melt is. 

However, the magmatic rocks which are widely used in 
fiber manufacture are more complicated multicomponent 
systems which include at least nine oxides, each of which 
can affect the viscosity to some extent. In addition, as numer¬ 
ous experimental studies of raw materials from different de¬ 
posits have shown, the effect of any particular oxide on the 
melt viscosity is determined not only by its nature but also by 
its content. Silica Si0 2 , alumina A1 2 0 3 , magnesium oxide 
MgO, and trivalent iron Fe 2 0 3 are viscosity-increasing oxi¬ 
des. Alkali metal oxides K 2 0 and Na 2 0 and divalent iron 
FeO lower the viscosity; the effect of others, including CaO, 
is quite complex and not unique: up to certain content they 
decrease the viscosity and above this value the viscosity of 
the melt increases or vice versa. 
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Repeated attempts to predict the viscosity of melts from 
the chemical composition of the rocks have been made. 
Using of mathematical modeling the authors of [1] obtained 
for two sets of raw materials with different chemical compo¬ 
sition regression equations relating the values of the loga¬ 
rithm of the viscosity at temperature 1400°C with the content 
of the main oxides in the rock: 

log t| 140 o o c = 0.658 + 0.003Si0 2 - 0.0057Ti0 2 + 
0.0025A1 2 0 3 - 0.0078Fe 2 0 3 - 0.0039FeO + 
0.0028MgO - 0.0027CaO - 
0.0003Na 2 0 + 0.0036K 2 0; (1) 

log t| 1400 o C = 1.613 + 0.0023Si0 2 - 0.0038Ti0 2 + 
0.0023A1 2 0 3 - 0.0094Fe 2 0 3 - 
0.0046FeO + 0.003 lMgO - 
0.0036CaO - 0.0014Na 2 0 + 0.00226K 2 0; (2) 

where Si0 2 , Ti0 2 , and others denote the oxide content, molar 
fraction x 10 3 . 

The equation (1) was derived for rock melts for which 
log r| varies from 0.5 to 3.6, while in the second equation this 
quantity is in the range 1.4 - 2.5. The coefficients of the fac¬ 
tors in the equations are different; correspondingly, the vis¬ 
cosity values obtained from Eq. (1) and (2) are also different. 
In addition, for the chosen composition the individual oxides 
MgO and K 2 0 increase the viscosity according to the first 
equation and decrease it according to the second equation. 

Similar relations have also been obtained for other sets of 
raw materials at different temperatures. In addition, the com¬ 
putational accuracy for log r| at 1450 and 1250°C is much 
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lower than at 1300 and 1350°C, while the change in the coef¬ 
ficients in these equations from one temperature to another 
has no logical explanation. On this basis the equations pre¬ 
sented cannot be regarded as correct and the viscosity values 
calculated using them as sufficiently accurate predictions of 
its real values. 

The author of [2] proposed a relation, which includes the 
acidity modulus M a , that he derived from experimental re¬ 
sults obtained for batch melts with different A1 2 0 3 , MgO, 
and CaO contents: 


Another method with a similar drawback for calculating 
the viscosity of magmatic melts from the chemical composi¬ 
tion of the rock is Persikov’s semi-empirical method [8]: 


log 11 = log 


(a-bK)x. 10 3 

4.576r-3.5+a(/? tot -p HiG )’ 


where r| is the viscosity, P; a and b are empirical regression 
coefficients; K is a structural-chemical parameter reflecting 
the ratio of the bridge 0° and non-bridge 0“ oxygen: 


log r| = 0.785M a - 0.00356^ + 5.11. 

Here M a = (Si0 2 + 2A1 2 0 3 )/(2Fe 2 0 3 + FeO + CaO + 
MgO + 2Na 2 0 + 2K 2 0 + MnO + Ti0 2 + ...), and t is the tem¬ 
perature, °C. 

The viscosity values computed at different temperatures 
and presented in [2] are largely at variance with the experi¬ 
mental data: the relative error fluctuates from 1.5 to 65%. 

A structural-thermophysical theory has been developed 
for magmatic melts [3]. According to this theory the visco¬ 
sity of melts can be calculated at their liquidus temperature 
without using the Frenkel-Andrade model theory of the vis¬ 
cosity [4]. The authors proposed using the thermophysical 
Fourier (Fo) and Biot (Bi) numbers and the Tikhonov num¬ 
ber (Ti) to evaluate the temperature dependence of the vis¬ 
cosity [5]. Assuming the thermal diffusivity and heat emis¬ 
sion for magmatic melts to be close to one another, they ob¬ 
tained the Einstein relation for the diffusion law in a different 
form [6] that made it possible to calculate the temperature 
variation of the melt viscosity as a function of the dimen¬ 
sionless numbers listed above: 

RT RT 

p - -== -, 

6nA Bi Fo VFo 6nA Ti Fo 

where p is the viscosity, P; R is a constant; T is the tempera¬ 
ture, K; A is the electronic polarization rate of non-bridge 
oxygen; and, Ti is the Tikhonov number (the homochronicity 
number for semi-infinite bodies). 

A nomogram representing the relation between the de¬ 
gree of polarization of non-bridge oxygen in melt and the 
Tikhonov number is used to calculate the numbers Bi and Ti 
for different values of the Fourier number [7]. 

Making definite suppositions and assumptions about the 
degree of polarization of non-bridge oxygen and the degree 
of depolymerization of the melt, the authors of [3] propose 
the following equation for calculating the viscosity at the 
liquidus temperature: 


log r\ = log 


RT 

671 4 Bi Fo VFo 


The viscosity values obtained with this equation likewise 
agree poorly with the experimental data. 


K = 


0“ 2(0-2M) 


O 


o 


M 


x 100: 


O and M are the total numbers of gram-ions of oxygen and 
network-formers (Si 4+ , Al 3+ , Fe +3 , P 5+ ) in the melt, respec¬ 
tively; T is the temperature, K; p is the pressure, MPa; a is 
the piezo-viscosity coefficient, MPa -i ; and, a= - 5.02 x 
10 ~ 4 MPa -1 for “dry” melts and a = -1.2x 10“'MPa -1 
for incompletely water-saturated melts. 

It should be noted that for calculations the author uses 
the concentrations of only three melt-forming oxides (Si0 2 , 
A1 2 0 3 , Fe 2 0 3 ), while their relative effects on the viscosity 
are neglected. 

In summary, no computational method studied above 
gives reliable values for the viscosity of fused glass and the 
methods are quite complicated to use in production practice. 
In this connection, using a multifactor regression analysis 
program the present authors attempted to obtain the best fit 
representing the dependence of the viscosity of melts on the 
chemical composition of the raw material. An array of expe¬ 
rimental data on the viscosity of basaltic melts of rocks from 
27 deposits with different chemical composition in the tem¬ 
perature range 1200 - 1450°C was used for modeling. 

The chemical composition of the rocks was determined 
by x-ray fluorescence using an SRM-25 computer-controlled 
multichannel x-ray spectrometer. Sample preparation was 
based on melting rock with active flux. State standard sam¬ 
ples of natural rocks were used to check the validity of the 
analysis. The determination errors of the concentrations of 
different oxides in the control rocks with confidence proba¬ 
bility 0.95 are presented in Table 1. The data are presented 
according to the OST 41-08212-82 second category of accu¬ 
racy. 

The viscosity at temperatures from 1200 to 1450°C was 
determined with a RVTS-K90R4 rotary viscosimeter. Its 
principle of operation reduces to measuring the time over 
which a coiled spring placed between a drive shaft and an in¬ 
ner cylinder, deforming the melt located between a stationary 
outer cylinder and the inner cylinder rotating with constant 
speed, turns through an angle. The absolute error of the mea¬ 
surement was ±0.5 Pa • sec. 

The temperature t , acidity modulus M a , and the mass 
content of oxides in basaltic rocks served as the initial pa¬ 
rameters Z t for obtaining an approximate relation. 
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TABLE 1. Determination Error for the Oxide Concentrations in the 


Control Rock 


Oxide 

Determination accuracy for standard samples 

Content, wt.% 

Absolute error, % 

Si0 2 

60.60 

±0.20 

A1 2 0 3 

18.20 

±0.20 

Ti0 2 

0.95 

±0.07 

MnO 

0.048 

±0.005 

Fe 2 0 3 + FeO 

7.20 

±0.15 

MgO 

2.20 

±0.10 

Na 2 0 

2.30 

±0.10 

k 2 o 

3.45 

±0.15 

p 2 o 5 

0.18 

±0.01 

Other 

4.70 

±0.20 


Note. Other — calcination losses. 


A preliminary analysis with 12 parameters showed that 
Si0 2 , A1 2 0 3 , CaO, and (FeO + Fe 2 0 3 ) have the strongest ef¬ 
fect on the melt viscosity. The effect of all other oxides is so 
weak as to be negligible. 

Assuming that the indicated parameters interfere with 
one another weakly, the approximate relation has the gener¬ 
ally accepted form 


hcalc > Pa • Sec 



ricaic, Pa • sec 



^ =f(Z l ,Z 2 ,Z 3 ,Z 4 ,Z 5 ,Z 6 ) = 

ydy yd 2 ya^ y a A y a 5 y a 6 
w 0 % ^2 3 "4 3 ^6 ’ 



Fig. 1. Relation between the computed and measured values of the 
melt viscosity at temperatures 1350 (a) and 1400°C ( b ). 


where Z 1 is the Si0 2 content, wt.%, Z 7 the A1 2 0 3 content, 
wt.%, Z 3 the CaO content, wt.%, Z 4 the (FeO + Fe 2 0 3 ) con¬ 
tent, wt.%; Z 5 is M a ; and, Z 6 is the temperature, °C. 

The problem was to find the seven unknown coefficients 
a 0 ... a 6 in the relation (3). 

The solution was obtained in two steps. 

First, the least-squares method was used to find the coef¬ 
ficients to a first approximation (a 0 , a { , ..., a 6 ), after which 
using direct-search technology the values were refined tak¬ 
ing account of the relative effect of the parameters on the vis¬ 
cosity. 

At the first step, denoting the experimental value of the 
viscosity of basaltic rocks as r| • (J = 1,2, ..., 27) and the vis¬ 
cosity values computed with the relation (3), into which the 
experimental data on the parameters Z t (/ = 1, 2, ..., 6) are 
substituted, for theyth deposit as 

Tl(fl 0 Z l) Z 2 j z y Z Tj Z 5j Z Tj ) = *1 (Zy X 


The sum of the squares of these differences is 

F ^Z F J = S[Uy -n(Z,)] 2 . 

7=1 


Next the minimum of the function F s was found with re¬ 
spect to each parameter a t . This gave the following system 
of seven equations with seven unknowns a 0 , a h ..., a 6 : 


d \F± 

d a 0 

d 

d a { 


= 0 
= 0 



d F± 

da 6 



the squared differences between the experiment value r|. of 
the viscosity and the computed value r|(Z 7 ) were written as 

IJ s 

Fj = [r( / - r|(Z..)] 2 . 


The coefficients a i to a first approximation were deter¬ 
mined by solving the system of equations (4) numerically. 

At the second step refined values were found for the co¬ 
efficients a i using a computer program to perform a direct 
search. 
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For this we wrote down the equation 

r\ = a 0 Z * 1 Z 2 2 Z 3 3 Z " 4 Z a 5 5 Z a 6 6 . (5) 

Each coefficient a i was varied over the range 
0.8a; < a i < 1.2 a i for all 27 rocks and the following minimum 
was found by the direct calculation using the relation (5): 


It is especially important to use a computational method 
when using in fiber manufacture rocks whose melts reach the 
required viscosity at temperatures above 1450°C and techni¬ 
cal difficulties appear in the experimental determination of 
the viscosity. 
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The regression analysis established that there a real rela¬ 
tion exists between the chemical composition of the rocks 
and the melt viscosity. This relation is described by a 
multifactor regression equation reflecting the individual and 
total M a contribution of the melt-forming oxides: 
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A positive sign of the exponent of oxide mass quantities 
signifies that for this set these exponents increase the visco¬ 
sity; a negative sign signifies a decrease. 

The equation obtained makes it possible to predict the 
melt viscosity to adequate accuracy at a prescribed tempera¬ 
ture from the chemical composition of the raw material. This 
is confirmed by the data presented in Fig. 1 for rocks which 
were not included in the regression analysis. 
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